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I. INTRODUCTION 

Over the last few years the Gutzwiller variational the- 
ory has developed into a useful tool for correlated multi- 
band systemsii^ii^ii^'^'i'^'^. Gutzwiller introduced his 
wave function to study ferromagnetism in the one-band 
Hubbard modeli^. For the evaluation of expectation val- 
ues he used a classical counting scheme, the so-called 
'Gutzwiller-approximation'; see also Re f. [ill for the dis- 
cussion of its physical content and Ref. [l2l for its math- 
ematical formulation. It was found later that this some- 
what ad hoc approximation is equivalent to an exact 
evaluation of expectation values in the limit of infinite 
spatial dimensions or lattice coordination number— 
An evaluation of Gutzwiller wave-functions in this limit 
will be denoted the 'Gutzwiller-theory' throughout this 
work. The limit of infinite spatial dimensions is also the 
central assumption in the dynamical mean-field theory 
(DMFT)i^. The solution of the DMFT equations is quite 
challenging even for the one-band Hubbard model. In 
contrast, the corresponding Gutzwiller variational space 
is just one-dimensional and its minimisation is a triv- 
ial numerical task. Of course, the study of the elec- 
tronic properties of real materials requires the treatment 
of multi-band Hubbard models. This task was accom- 
plished some years ago in Refs. [l]I3[l]. 

An alternative scheme to derive the Gutzwiller energy 
functional for a single band is based on the slave-boson 
mean-field theory (SBMFT) of Kotliar and Rucken- 
steiriF^. A generalization of this approach which repro- 
duces the results of the multi-band Gutzwiller theory 
is straightforward for systems with only density-density 
interactiortiiiiiii^. It was only recently, however, that 
Lechermann et al.— succeeded to develop a SBMFT 
scheme that allows to investigate systems with general 
multi-band interactions. The authors interpret their ap- 
proach as a generalization of the Gutzwiller theory. In 
fact, as we will show in this work, both theories are com- 
pletely equivalent. 

Our article is structured as follows. In Sect. |TT] we 
introduce the multi-band Hamiltonian and the general 
class of Gutzwiller wave- functions. The equivalence of 
the SBMFT, as derived in Ref. [li|, and the Gutzwiller 
theory is demonstrated in Sect. IIIII Finally, we critically 
discuss the limitations of both methods in Sect. IIVI 



II. HUBBARD MODELS AND GUTZWILLER 
WAVE-FUNCTIONS 

We investigate multi-band Hubbard models, described 
by the general class of Hamiltonians 

^= E </'Cs-,.'+E^i-.^ = ^o+^i-. (1) 

Here, the first term describes the hopping of electrons be- 
tween spin-orbital states a, a' on lattice sites respec- 
tively. The Hamiltonian H\oc,i contains all local terms, 
i.e., the two-particle Coulomb interactions and the or- 
bital onsite-energies. For any lattice site i one introduces 
the Fock-states \I)i, in which certain sets of spin-orbital 
states a are occupied^ii. These states form a basis of the 
local atomic Hilbert space and can be used to write any 
other local multiplet state as 

|rK = E^/Tl^)- (2) 

/ 

The most general Ansatz for a multi-band Gutzwiller 
wave-function has the form 

i 

where is a normalized single-particle product state 
and the local Gutzwiller correlator is defined as 

A = E4T'|r>..(r'|. (4) 

r,r' 

Here, the states \T)i can be an arbitrary atomic basis 
and the numbers ApV' are variational parameters. In 
our first work, Ref. assumed that the multiplets 

\T)i are the eigenstates of H\oc,i and App, ^ (5r.r'- The 
more general Ansatz (|4]) was first evaluated in Ref. Q for 
Hermitian operators Pg,j:- The non-Hermitian case has 
been studied in Ref. 8]. In the following, we drop the 
site index when we deal with purely local quantities. 
In general, the uncorrelated local density matrix 

Cl,, ^ {cU..ho (5) 
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is not diagonal. It is then useful to introduce a second 
orbital basis, defined by the operators h':^\ which, by 
construction, have a diagonal local density matrix, 

H^' ^ E "x^^^^ ' {KK')'^o = Ky< ■ (6) 

cr 

Within the Gutzwiller theory one usually works in the 
new '/i-representation' because all formulae have a much 
simpler form than in the original 'c-representation'. How- 
ever, in order to show the equivalence of the Gutzwiller 
theory with the slave-boson results in section HTTl we have 
to work with both representations simultaneously. The 
operators hl^ define Fock states \H) which can also be 
used to write the multiplet states ^ as 

\T)^J2TH,r\H), (7) 

H 

where the coefficients Tn.r and T/.r are related through 
Thx = E ^HjTi^r , ^Hj ^ {H\I) . (8) 

III. COMPARISON 

The calculation of expectation values for Gutzwiller 
wave- functions is a straightforward task, once the ba- 
sic structure of diagrams in infinite dimensions is under- 
stoodii^i^. In contrast, the derivation of the generalized 
slave-boson mean-field theory in Ref. ,19] requires a num- 
ber of subtle ideas. Furthermore, there is not a clear 
correspondence of all the mathematical objects that ap- 
pear in both approaches. Therefore, we are not going 
to compare any particular steps of the two quite differ- 
ent derivations, but focus on the final energy functional 
at zero temperature in order to show the equivalence of 
both approaches. For all details of the derivations we 
refer the reader to Refs. jL>2.,4.19.]. 

A. Local energy 

In infinite dimensions, the expectation value of the lo- 
cal Haniiltonian H\oc,i in the Gutzwiller wave function 
reads 

(Aoc,j)*G = E '^r2,ri'^r3,r4-^rT,r3"^ri,r4 ' (9) 
ri...r4 

where 

i?^-P3 ^ (r2|i/ioc..|r3) , (10) 
m°^,P^ ^ ((|ri)(r4|))*„ . (11) 

The expectation value (jlip can be written as 

'^riT4 = E ^Hi,ri^H4,r4"^ffi,-ff4 • (12) 

Hi ,//4 

= J2^H,rXHX,< (13) 

H 



because, for a diagonal local density-matrix in the h- 
representation, one readily finds 

'rriRM' = SH.H'm^H ' (14) 

^ n < n (i-o- (15) 

7(occ.) 7'{unocc.) 

In order to make contact with the results in Ref. [l^ we 
need to bring (jl3p in the slightly more complicated form 

"^ri,r4 = E '^h,Ti^'hj\P^'^h',tS^H'j\P^' 

H,H',I 

(16) 

which is equivalent to because of the completeness 
relation 

Y,^H'jn*Hj = 5H,W ■ (17) 
/ 

We now introduce the new variational parameters 

tprj = ><TX'Thx'^Hj\fm^ (18) 

T',H 

which allow us to write the expectation value ^ as 

(^loc.)*o = E E ^tjVt'jE'^x' (19) 
r,r' / 

This equation has exactly the same form as equation 
(47) in Ref. [l^, after the slave-boson operators 0r./ 
have been replaced by their mean-field expectation val- 
ues, (j)T,i ^ ipr,i- 

B. Local constraints 

The variational parameters need to obey certain con- 
straints which naturally arise in the evaluation in infinite 
dimensions. These are 

(Ptp)*„ ^ 1 , (20) 
{cic,,P^PU„ = (45.')*o ■ (21) 

Note that moving the operator P^p relative to cj. or c^, 
in ([2T|l would not alter the whole set of constraints. A 
set of constraints equivalent to (|2ip is obtained when we 
use the operators hj'^ , 

{h\h^,P^P)^„ - {hlh^,h„ . (22) 
The constraint (|^n|) can be written as 

E '^r,riAr,r2™ri,r2 = 1 > (23) 
r,ri,r2 

which, by use of eqs. and is found to be equiv- 
alent to 

J2'Pr,i^rj = '^- (24) 
r,/ 
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This is equation (28) in Ref. [19| at mean- field level. For 
the constraints (l22ll it follows that 



(25) 



r,ri H\,H2 



where, due to eqs. and ([TS]), the expectation value 
on the l.h.s. can be written as 



Then, the identity 



(26) 



{H2\h\h^.\H^) = ^H.M{h\h\h^>\h)n*H,j, (27) 



transforms eq. (|25p into the form 

E E ^T,iVT,r {I'lhlhy \I) = {h\k,,h, . (28) 



r /,/' 



These equations can be transformed to the c-representa- 
tion which leads to 

EE^rj'^r,7(/|c^c,,|/') = (ct . (29) 
r /,/' 

Equation ([29|) is equivalent to equation (29) in Ref. [l^ 
at mean- field level. 



C. Hopping renormalization 

Finally, we investigate the expectation value of the 
electron transfer operators in the Hamiltonian ([T]). In 
infinite dimensions one finds that such an expectation 
value has the form (i ^ j) 

or, alternatively, in the /i-representation 

C^lAy^o = E (\U;^^,^>*o ■ (31) 

The local renormalization-matrix is most easily calcu- 
lated in the /i-representation^ where it has the rather 
simple form 



(32) 



The matrix in the c-representation can then be de- 
rived from q"^ by the transformation 



(33) 



With eqs. ([3]) and (|4]), the matrix q^ reads explicitly 

97' = Jr E ^r2.r,Ar3.r.(r2|/.;|r3) (34) 
7' ri...r4 

x((|r,)(r.|v))^^ ' 

where the expectation value in p4p can be written as 
(|ri)(r4|v))^^ = E TH,,r^T^,,rAH,\hy\H,) 

Hi ,H4 



x,/m%ml^^^-L^. (35) 



With eqs. p?]) and (fTSl) we can rewrite p4|) as 



,0(/_„0)EE^r./^r',/' 



x(r|/.t|r')(/'IV|/). (36) 

The transformation (j33p from the /i-representation to 
the c-representation is not as straightforward as the cor- 
responding transformation from eq. to eq. (|29p . 
Whereas the transformation with respect to the lower 
index 7 is still simple 



(r|ct|r') = E(r|/^;ir> 



(37) 



for the upper index 7' we need to take into account the 

factor ^l/{ny{l — riy)) in ((36|) which also depends on 

7'. For this purpose, we introduce the hole density- 
matrix D with the elements 

D^,, EE {c,,ciU, , (38) 

in addition to the density-matrix C already defined in 
([5]) . Then the transformation ([55]) for the upper index 7' 
can be carried out along the lines 



E"7'-' 



E 



EE E ((C°^°)-^ 



(39) 



CT.cr' 



Here, we used the notation 



7' V4(l""7') 



(40) 



Witheq. ^ and with (/'|c^|/) = (/|ct|/') we can finally 
write the renormalization-matrix in the c-representation 
as 



7,7' 



EE ^r,7^r',7'(r|4|r') 
r,r' I J' 

xE((<5°^°)-^). (41) 
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This expression matches equation (37) in Ref. [T^ at 
mean field level, apart from the fact that there the con- 
straints have been used to write the matrices C^D° as a 
function of the fields (^p j ; see equations (35) and (36) in 
Ref. pjj]. However, as long as the constraints are fulfilled, 
this makes no diff'erence because it does not change the 
variational energy functional. 

IV. DISCUSSION 

In this work we showed that the multi-band slave- 
boson mean-field theory of Lcchcrmann et al. repro- 
duces correctly the energy functional of the multi-band 
Gutzwiller theory developed earlier. As a byproduct we 
were able to show how the different objects that appear 
in both approaches are related. This will turn out to be 
important for the comparison of future numerical results. 

We believe that there are good reasons to prefer the 
derivation based on Gutzwiller wave-functions over the 
slave-boson mean- field theory. First, these wave func- 
tions are well defined and they are evaluated exactly 
in the unambiguous limit of infinite spatial dimensions 
(D — > oo). Therefore, e.g., the inclusion of superconduct- 
ing pair-correlations was straightforward^'*'. In contrast, 
the slave-boson mean-field derivation is uncontrolled and 
quite adjustable in its outcome. Of all the different equa- 
tions that one may derive within such an approach, the 
'right ones' are usually identified by some sophisticated 
guess. This guess, not surprisingly, always turns out to 
be equivalent to the Gutzwiller theory. This very equiv- 
alence is, by far, the most convincing argument for the 



credibility of the SBMFT results. 

The Gutzwiller theory can also be used to calcu- 
late quasi-particle excitations within a Fermi-liquid ap- 
proact4i^ as well as spin-wave excitations^. The 
quasi-particle bands in the Gutzwiller theory coincide 
with those derived in the SBMFT. Therefore the zero- 
temperature spectral properties are equivalent in both 
approaches. 

The ground-state energy functional also provides the 
Landau parameters for the description of thermodynamic 
properties'^. Therefore, both approaches are equivalent 
in the Fermi-liquid regime when the temperatures T is 
much smaller than the Fermi temperature Tp. Although 
the SBMFT equations can also be solved for T « Tp or 
even for T ^ Tp, the approximation breaks down in this 
temperature regime^.-^. 

In principle, the SBMFT could be improved by com- 
puting fluctuations around the saddle point. To the best 
of our knowledge, however, this promise, although often 
made, has never materialized in any convincing improve- 
ment of the results, not even for the one-band Hubbard 
modeP'^. In contrast, for Gutzwiller wave- functions it is 
possible to calculate systematically 1/D corrections for 
all physical quantitiesi^^. Such calculations allow to es- 
timate the accuracy of the results in infinite dimensions 
and to improve them, if necessary. 

It should be kept in mind that the Gutzwiller theory 
is based on rather simple variational many-body wave 
functions which could be improved in many directions. 
Despite its limitations, however, the Gutzwiller theory 
appears to provide a suitable description of the quasi- 
particle bands in ferromagnetic Nickel^i^. 
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